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^SJ ' We present the photometric calibration of the Swift UhraViolet/Optical Telescope 

(UVOT) which includes: optimum photometric and background apertures, effective 
area curves, colour transformations, conversion factors for count rates to flux, and the 
photometric zero points (which are accurate to better than 4 per cent) for each of 
f-*) I the seven UVOT broadband filters. The calibration was performed with observations 

• • . of standard stars and standard star fields that represent a wide range of spectral 

' star types. The calibration results include the position dependent uniformity, and 

v> I instrument response over the I6OO-8OOOA operational range. Because the UVOT is 

j_j ■ a photon counting instrument, we also discuss the effect of coincidence loss on the 

' calibration results. We provide practical guidelines for using the calibration in UVOT 

data analysis. The results presented here supersede previous calibration results. 
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1 INTRODUCTION 



The Ultraviolet Optical Telescope (UVOT; iRoming etHI 
I2OO5I ) is a modified Ritchey-Chretien telescope with a 
17 X 17 arcmin field-of-view (FOV) operating in the 1600- 
8OOOA range. Like many optical telescopes, the UVOT uses 
broadband filters in the ultraviolet (UV) and optical to ob- 
tain colour information. However, unlike most optical tele- 
scopes, the UVOT has a photon counting detector that gath- 
ers data in a similar way to an X-ray detector. 

The UVOT is o ne of three telescope s flying on board 
the Swift spacecraft (jCehrels et al. |[2004l ) and is co-aligned 
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with the 15 - 150 k eV Burst Alert Telescope (BAT; 
iBarthelmv et al~ II2OO5I) and the .2 - 10 keV X-Ray Tele- 
scope (XRT; iBurrows et aTlbOOSl l . The primary goal of the 
Swtft mission is to detect and characterize gamma-ray bursts 
(GRBs) and their afterglows. The design of the UVOT is 
well suited to this goal. In addition, UVOT's UV response 
in particular makes it a valuable instrument for other types 
of observations. 

There are currently four near-UV imaging telescopes 
operating in s pace: HST - Wide-F ield Planetary Cam- 
era 2 (WFP C2: iBurrows et al~lll994h . VMM-Optical Mon- 
itor (OM; iMason et al. I I2001I'). Galaxy Evolution Ex- 
plorer fCALEX: [siegmund et all 12004 iMiUiard et "aTl 
I2OOII : iBianchilbOOd ). and Swift - UVOT. The HST has the 
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highest spatial resolution and sensitivity but the smallest 
FOV. GALEX has the largest FOV and covers a larger UV 
wavelength range than OM or UVOT but it has very broad- 
band filters and the lowest spatial resolution of the four tele- 
scopes. The UVOT and OM are very similar in design, both 
covering approximately the same wavelength range and con- 
taining similar broadband filters in both the UV and optical 
regimes. Both UVOT and OM include three UV fihers that 
together cover approximately the same band-pass as the long 
wavelength GALEX fiher. The UVOT has a shghtly broader 
point spread function than OM but a higher sensitivity, by 
a factor of ~ 10 at the bluest wavelengths. With respect to 
the four telescopes, UVOT, being mounted on Swift, is able 
to respond most quickly, making it the best instrument for 
observing transient phenomena. 

In this paper, we describe the in-orbit photometric cal- 
ibration of the UVOT and define the UVOT photometric 
system. We begin by describing the UVOT technical details 
in Section [21 the calibration strategy in Section O followed 
by the observational measurements and photometry method 
in Sections|3HSl We review the calibration of the coincidence 
loss and position-dependent uniformity in Sections [7] and |8l 
respectively, and we provide the in-orbit effective area cal- 
culations, photometric zero points, colour transformations, 
and flux conversion factors in Sections l9l 412l We test the 
calibration in Section [13] and give guidelines for making use 
of the calibration for UVOT data analysis in Section [HI Fi- 
nally we summarise the calibration and discuss what is left 
to be done in Section [TSl 

Our preliminary UV OT calibration, d e scribe d in 
[Breeveld et al. I l|2005l ) and [ivanushkina et al. I (|2005l ). was 
released soon after launch in the HEASARC calibration 
database (versio n 20050805 of the Swift/UVOTA CalD^. 
[Li et al.l (|2006h later gave an independent calibration for 
the optical filters using data from early in the mission. The 
calibration analysis described in this paper is the most de- 
tailed and comprehensive to date and thus supersedes our 
previous work. 



2 UVOT TECHNICAL DETAILS 

A technical descrip tion of the UVOT is given by 
[Roming et al. [ l|2005h : we summarise the relevant properties 
here. The UVOT is of a modified Ritchey-Chretien design, 
with a 30cm primary mirror and an /-ratio of //12.7 after 
the secondary. A 45 degree mirror directs photons into one of 
two detector units, one of which is kept in cold redundancy. 
Each detector unit consists of a filter wheel and micro- 
channel plate intensifi ed CCD (MIC: lKawakami et al. [[19941 : 
[Fordham et al. |[l992l ). The filter wheels contain 11 slots 
which house 3 optical filters, 3 UV filters, a clear white light 
filter (which transmits throughout the UVOT wavelength 
range of I6OO-8OOOA), a magnifier, a low-resolution optical 
grism, a low-resolution UV grism, and a blocked filter. The 
characteristics of the optical and UV filters can be found in 
Table [1] 

The MIC detector has an entrance window that is 
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Table 1. Swift/VYOT Filter Characteristics. The central wave- 
length is the midpoint between the wavelengths at half maximum 



Fiher 


Central Wavelength 

(A) 


FWHM 

(A) 


V 


5468 


769 


b 


4392 


975 


u 


3465 


785 


uvwl 


2600 


693 


uvm2 


2246 


498 


uvw2 


1928 


657 



slightly figured optically to flatten the image plane of the 
telescope. An S20 photocathode is deposited on the inside 
of the window and is optimized for the UV and blue wave- 
lengths. The photocathode converts an incoming photon 
into an electron signal, which is then amplified by a factor of 
a million by a photomultiplier stage. Thus every incoming 
photon results in a cloud of electrons at the back of the pho- 
tomultiplier, and these are converted back into photons by 
a phosphor screen. The resulting photon splash is recorded 
on a fast scan CCD, which has an active area of 256 x 256 
pixels. Each photon splash extends over several CCD pixels, 
and this allows the centroid of the splash to be determined 
in real time by fast onboard electronics to a fraction of a 
CCD pixel (1/8*' ' of a CCD pixel in the case of the UVOT 
MIC detectors; [Michel et aLl[l997[ '). The resuhing image 
format is thus 2048 x 2048 pixels, covering a field of view of 
17 X 17 arcmin. This provides a spatial sampling of 0.5 arc- 
sec. There is a fibre taper between the phosphor screen and 
the CCD, which compensates for the larger physical area of 
the photomultiplier stage compared to the CCD. Through- 
out this paper, the word pixel refers to 0.5 x 0.5 arcsec image 
pixels, unless explicitly stated otherwise. 

The advantages of the MIC detector over a traditional 
'bare' CCD are that (a) it is photon counting, (b) it op- 
erates at 'room temperature' without the need for cooling 
and (c) it is insensitive to cosmic ray and charged particle 
hits in the CCD. Because it is a photon counting detector, 
the UVOT suffers from coincidence losses at high photon 
rates, when two or more photons arrive at a similar loca- 
tion on the detector w ithin the same CCD read out interval 
([Fordham et al. [[2OO OY It is equivalent to 'pile-up' in an X- 
ray CCD detector. The magnitude of the effect depends on 
the read-out rate of the CCD, which is once every 11.0329 
ms when the full 256 x 256 CCD pixels are used, and faster 
when only part of the field is read out. In addition to coin- 
cidence loss the count rates also have to be corrected for the 
dead time while the charge is transferred out of the CCD, 
which amounts to 1.6 per cent of the full-frame read-out in- 
terval. The loss due to dead time is a constant factor for a 
given read-out rate, whereas coincidence loss also depends 
on the incoming photon rate. 

For the full CCD, with a read-out time of ~ 11 ms, 
coincidence losses start to be significant at ^ 10 counts s~^ 
and a correction should be applied to the recorded signal. 
Beyond ~ 90 counts s~^, equivalent to 1 count per frame, 
correcting for coincidence losses becomes increasingly uncer- 
tain. To extend the range of the detector to brighter sources, 
the read-out rate can be increased, as noted above, by using 
only a subset of the CCD pixels. We use two pre-defined 
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hardware windows, which are 120 x 120 CCD pixels (8x8 
arcmin) and 75 x 75 CCD pixels (5x5 arcmin) on a side, cen- 
tred on the observatory boresight. These reduce the frame 
times to 5.417ms and 3.600ms, respectively, and also reduce 
the dead time by a small amount. 

The onboard algorithm that centroids the photon splash 
recorded on the MIC CCD uses a simple look-up table for 
speed. The intrinsic imperfections in this process mean that 
the sub-prxels within each physical CCD pixel do not all 
have the same effective physical area. The signature of this 
in the resulting image is an apparent 8x8 pixel fixed-pattern 
effect, sometimes referred to as 'mod-8 noise' , though it is 
not s trictly noise since photons are conserved (|Michel et al. I 
1 19971 ). A calibration of this is derived onboard based on illu- 
minating the field with an LED within the instrument, and 
applied as a correction to the centroiding algorithm. Nev- 
ertheless small gain variations over the face of the detector 
mean that a low-level residual fixed-pattern effect is still 
seen which varies in magnitude depending on position on 
the detector. If desired, the effect can be largely removed on 
the ground by using the tool uvotmodmap, which is released 
as part of the HEAsoft Swift ftools software packagfl 

The UVOT operates in two data-taking modes: event 
and image. Event mode preserves the temporal and posi- 
tional information of each photon. In image mode the pho- 
tons are accumulated into an image in the instrument's on- 
board memory before being telemetered to the ground, thus 
providing positional but not temporal information on the 
photons. 

Because the spacecraft can drift slightly during an ob- 
servation, the photon positions are shifted before they are 
added to memory. Occasionally, due to large spurious drifts 
in the spacecraft attitude information, individual frames are 
not added to the image. This time lost, when the onboard 
shift-and-add algorithm unnecessarily tosses events off the 
image, is known as TOSSLOSS. TOSSLOSS occurred fairly 
often until a software fix in September 2005, but scarcely 
ever since then. 



3 CALIBRATION STRATEGY 

In the following sections we describe the different elements 
of the UVOT calibration in detail, but first we will outline 
the overall approach that we have taken. As a starting point, 
we used our ground-based measurements of the various com- 
ponents in the UVOT optical path, to produce an idealised 
or reference effective area curve for each of the UVOT fil- 
ters. From these effective area curves, we predicted the count 
rates for a number of photometric standard stars for each 
of the UVOT filters. We then compared these predictions 
with in-orbit measurements of the standard stars, to pro- 
duce an in-orbit correction curve to the overall instrument 
response as a function of wavelength. Our predicted effec- 
tive area curves were multiplied by the in-orbit correction 
curve to obtain the in-orbit effective area curves. Once the 
in-orbit effective areas were established, zero points, colour 
transformations, fiux conversion factors etc. followed. 

^ HEAsoft software can be found at 

|http : / /heasarc . gsfc .nasa. gov /docs/ software /Iheasoft / 1 



Any changes in the overall filter transmission are thus 
incorporated into the correction curve. Changes in filter 
transmission shape would be more difficult to deal with, but 
would show up as a difference in zero point for sources of 
different colour. We do not find evidence for a change in the 
transmission shape for any of the filters. 



4 UVOT INSTRUMENT RESPONSE 

The instrument response of the UVOT is a product of: 

(i) the telescope primary mirror geometric collecting area 
of 659cm^, 

(ii) the mirror reflectivity (|Roming II2007I ). 

(iii) the filter transmission curves, 

(iv) the detector quantum efficiency (DQE), which is the 
overall sensitivity of the phot on counting syste m, including 
the photocathode sensitivity (|Kawakami Ill999h . 

The DQE and filter transmission curves (for all but the 
white filter, see below) as a function of wavelength were 
measured in the laboratory at a sub-unit level. The mirror 
reflectivity was also determined in the laboratory, by mea- 
suring planar witness samples that were coated alongside 
the mirrors. 

The measurements listed above, but not including the 
fllter transmission curves, were combined to produce an ide- 
alised, or 'reference' response for the telescope, which is 
shown as a dashed line in Figure [T] This is effectively the 
response expected if all the photons that pass through the 
telescope can be captured in the image. In practice, residual 
imperfections in the reflecting surfaces will scatter some pho- 
tons into broad wings in the telescope point spread function, 
and thus the actual throughput for photons in the image core 
will be less than ideal. 

The dotted line in Figure [T] is placed at the wavelength 
of I6OOA, which is the short wavelength limit of the instru- 
ment response determined by the detector window trans- 
mission. The long wavelength cut-off has been set at 8OOOA. 
Outside the limits of I6OOA and 8OOOA the detector quan- 
tum efficiency is very small; this range includes 99.98 per 
cent of the total instrument response. 

Figure [5] shows the transmission curves of the opti- 
cal a nd UV filters as measured in the laboratory (iRoming I 
l2003h . The extra peak at ~ 4750 A is part of the u fiher 
transmission. The error at each point on these transmis- 
sion curves is at most ±1 per cent (accounting for systemat- 
ics), but is probably ±0.4 per cent. The curves provided for 
uvw2 do not extend shortward of I8OOA due to the inability 
of the measuring device to provide data below this wave- 
length; therefore this filter curve had to be extrapolated to 
the I6OOA limit assuming a peak at 1930A, and a symmetri- 
cal response profile. There were no laboratory measurements 
of the white light filter transmission as a function of wave- 
length. Instead, the curve that we have adopted, and which 
is shown in Figure (2] is a model of the transmission, based 
on the design and substrate properties of the white filter. 

The predicted effective area curves for each of the filters 
were obtained by combining the instrument response curve 
(Figure [l]) with the individual filter transmission curves 
(Figure [U. 
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2000 4000 6000 BOOO 

Wovelength (A) 

Figure 1. The instrument response of the UVOT including de- 
tector window, cathode sensitivity, mirror reflectivity and tele- 
scope area. The dashed line represents the idealised instrument 
response (assuming all photons are collected in the image) while 
the solid line shows the actual in-orbit response measured in a 
5 arcsec radius circular aperture. The vertical dotted line marks 
the short wavelength cut-off of the instrument at I6OOA. 




4000 5000 6000 7000 8000 

Wavelength (A) 

Figure 2. Filter transmission curves for the UVOT filters as 
measured in the laboratory. The white filter transmission curve 
is given by the dashed line; the identities of the other filters are 
indicated on the plot. The extra peak at ~ 4750A is part of the 
u filter transmission. 



5 STANDARD STAR DATA 

Table [2] lists the standard stars and standard fields that 
were used for the calibration work presented in this paper. 
In total, two Oke standard stars, six white dwarf standards, 
ten Landolt standard stars, and four star fields have been 
used for the various calibration functions listed in Table [2] 
For the standard stars used to determine the in-orbit ef- 
fective areas, zero points, colour transformations and count 
rate to flux conversion factors, our calibration procedure re- 
quires spectra that can be folded through the instrument re- 
sponse and filter transmission curves. It has been a challenge 
to find spectrophotometric standard stars that are known to 
sufficient accuracy, have wide enough spectral coverage, are 
not variable, and are not too bright for the sensitive UVOT 




500C 
Wavelength (A) 



Figure 3. Comparison between the Johnson and UVOT 
optical responses. The Jo hnso n responses are taken from 
lAzusienis &: Straizvs I lli969l ) and iBuser I l[l97^ , as described in 
Section |5] 



detectors. Thus we have been limited in the UV to just 3 
trusted white dwarf standards. For the optical we have used 
photometric standards to increase the sample of sources. 

Representations of the Johnson UBV passbands are re- 
quired for several aspects or our calibration (e.g. normalisa- 
tion of standard star spectra as described below and colour 
transformations from UVOT to the Johnson system). How- 
ever, the true shapes of the J ohnson UBV response curves 
have long been controversial (|Bessel 1 12005| ). In this work, 
we have taken the following response curves to represent 
the Johnson UBV system. For the V and B responses we 
have use d Table 1, column s and "</>s" respectively, 

from Azusienis fc Straizvs I (Il969l ) . For th e John son U re- 
sponse, we have used Table 2 from iBuser I (Il978l ). Figure [3] 
shows the comparison between the Johnson U, B and V and 
UVOT normalised u, b and v responses. In this paper we 
use the convention of capital letters for the Johnson system 
magnitudes and lower case letters for the UVOT system 
magnitudes. 

For the two Oke standard stars, spectra were 
obtained from the ESO standard star archive 
(ftp:/ /ftp.eso.org/pub/stecf/standards/hststan ). These 
spectra cover the wavelength range 32OO-8OOOA and are 
therefore on ly used to calibrate the v and b filters. As 
described in IColina fc BohlirTI l|l994 ). the absolute photo- 
metric calibration of these spectra can be improved using 
precision ground-based photometry. These spectra were 
therefore folded through the Johnson V filter response and 
norma lised to the V magnitudes given by IColina fc Bohlin I 
(| 19941 ). 

Three of the white dwarf standards (WD1657-I-343, 
WD1026-K453, WD1121-hl45) are used for the effective 
areas and zero point calibration. Spectra were obtained 
from the HST MAST archive (http^/archive.stscLedu/} 
of WD1657+34:J3 and WD1026+45fl The spectrum of 
WD1026-I-453 had to be extrapolated longward of 5700A, 



3 MAST IDs: OSviOlOlO, OShlilOlO, 08hlil040, 08vl01030 
" MAST IDs: 08hf06040, 08hl06010, 08hi06020 



Photometric calibration of the swift ultraviolet/ optical telescope 5 



and so this spectrum was not used for calibration of the 
V filter. F or WD1121+145 we use the l UE s pectrum pre- 
sented bv lHolberg. Barstow fc B urleigh' (2003') and the op- 
tical spectrum from iMassev et al. (1988. ) . For these three 
sources, the Johnson magnitudes listed in Table [2] were ob- 
tained by folding their spectra through the Johnson UBV 

response curves. 

The Landolt standard stars (jLandolt Ill992f ) do not have 
accurate spectra available in the literature. Therefore, the 
Landolt standard sta rs were matche d to the stellar spectra 
from the catalogue of IPicklesI (|l998l ) which had the closest 
B-V and U-B colours. The best-matching Pickles spectra 
were then folded through the Johnson transmissi on curves 
and n ormalised to the U, B, or V magnitudes from lLandolt I 
l| 19921 ) for calibration of the UVOT u, b and v filters respec- 
tively. 



5.1 In-orbit measurements 

Table [2] gives observational details for all of the calibration 
sources and fields, including observation dates (Column 2), 
the observed filters (Column 3), and the photometric cali- 
bration product in which the source or field was used (Col- 
umn 4). Most of the observations were taken between launch 
and April 2005, during the calibration phase of the Swift 
mission. For more details about the observational informa- 
tion please re fer to the Calibration Database (CalDB) 
documentfO (|Poole |[2007al lbll3l. The original ground-based 
UBV photometry for the Landolt standar d stars was ob - 
tained using a 14 arcsec diameter aperture (|Landolt II 19921 ). 
which is not optimum for the UVOT calibration. We there- 
fore checked each of the UVOT images of Landolt standard 
stars to ensure that there were no other stars of sufficient 
brightness to offset the photometric calibration, within 7 
arcsec of the standard stars. 



6 PHOTOMETRY METHOD 

In this section we describe how and why we chose the pho- 
tometric aperture and background region to use for the cal- 
ibration analysis, and describe the analysis we undertook. 
We do not necessarily recommend following the same proce- 
dure for all UVOT data analysis and the reader is referred 
to Section (21 for more information on this. 

Much of this analysis was performed using in-house 
IDlS routines which could be readily tailored to the cali- 
bration tasks at hand, and have allowed us to test analysis 
methods and calibration products that have subsequently 
been implemented in the official UVOT HEAsoft software 
and HEASARC calibration database. 



6.1 Optimum photometric aperture 

One effect of coincidence loss (see Section [Tjl is that the 
shape of the point spread function (PSF) is slightly depen- 
dent on the source count rate. We wished to find an aper- 
ture enclosing the same percentage of the PSF regardless of 



^ |http : / /swift . gsfc . nasa. gov/docs/heas arc/ caldb /s wift / 1 
^ from ITT Visual Information Solutions 




5 1 

Radius (pixels) 

Figure 4. Coincidence-loss-corrected count rates within concen- 
tric radii, for Hz2 in the uvw2 filter, with a hardware window 
size of 2480 X 2480 (solid black line), 960 X 960 (dashed line), 
and 600 X 600 (dotted line). The error bars on the background 
subtracted data were determined using Poisson statistics 



count rate, if such an aperture existed. This was most eas- 
ily found by using the same source and changing the frame 
rate to mimic a change in count rate and then looking for 
a consistent result after correcting for coincidence loss (see 
Section [T] Equation [TJ. The use of different hardware win- 
dow sizes, and hence frame times, allows coincidence loss to 
be varied without changing any other parameters. There- 
fore, we define the optimum photometric aperture radius to 
be that which gives rise to the smallest variation in the en- 
closed energy fraction for a point source as a function of 
coincidence loss. This is the most convenient aperture to 
use for the calibration because it gives consistent results for 
photometry over a wide range of count rates. 

The optimum aperture was investigated using obser- 
vations of two isolated, bright stars GD128 and Hz2, for 
the default set of hardware windows with sides of length 
2048, 960 and 600 pixels. Hz2 was used for the optical and 
UV filters, but is too bright for the white filter, for which 
GD128 was used. The raw count rates were corrected using 
the theoretical coincidence loss equation only (see Section [71 
Equation [l}, which takes the frame time into consideration. 
Figure |4] shows an example of this for the uvw2 filter. The 
minimum of the RMS curve coincides with the optimum 
aperture radius. 

Table [3] shows the results from Hz2 (optical and UV 
filters) and GD128 (white filter) for the optimum aperture: 
the optical, UV and white filters have optimum aperture 
radn of 10.1 ± 1.3 pixels, 10.7 ± 1.5 pixels, and 9.8 ± 1.3 
pixels respectively. Since these are similar, and consistent 
within the uncertainties with a mean of 10.5 ± 1.2 pixels, 
we chose to use a radius of 10 pixels (5 arcsec) for all the 
photometric calibration analysis. 

The percentage of counts within the chosen 10 pixel 
aperture radius is stable for moderate count rates. Assum- 
ing that a 55 pixel radius aper ture represents 100 per cent of 
the PSF (see Section [6.31 and iBreeveld etal"! (|2007D ). then 
85.8±3.8 per cent of the PSF is contained within the 10 pixel 
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Table 2. Observational data used for the UVOT photometric calibration. The uses are given in column 4 where: OA is optimum photo- 
metric aperture (Section l6.H l. BA is background aperture (Section l6.2l fc l63)l . CL is coincidence loss empirical correction (Section^, PDU 
is position dependent uniformity (Section [Sjl, EA is in-orbit effective area curve (Section |9]l, ZP is photometric zero points (Section IIOI I. 
and CT is colour tran sformations (S ection | ll[. The las t three columns show the Johnson V, B and U magnitudes; for the Landolt sources 
from iLandolt ] lll992h . the Oke from lColina fc Bohlin I ll 19941 ) and white dwarfs calculated as described in Section [5] 



Source 


Date 


Filter 


Use 


Origin 


V 




B 




U 




SA95-42 


July 05 


V 






EA, ZP, CT 


Oke 


15, 


.606 


15, 


.391 


- 




g24-9 


July 05 


V, 


b 




EA, ZP, CF 


Oke 


15, 


.751 


16, 


.176 


- 




WD1657-I-343 


Feb- April 05, June 05 


V, 


b. 


u, uvwl, uvm2, uvw2, white 


EA, ZP 


White Dwarf 


16, 


.4 


16, 


.2 


15. 


,0 


WD1026-I-453 


July 05, Oct-Nov 05 


b. 


u. 


uvwl, uvm2, uvw2 


EA, ZP 


White Dwarf 


16, 


.1 


15, 


.9 


14. 


.8 


WD1121-I-145 


Feb-May 05 


uvwl, uvm2, uvw2, white 


BA, EA, ZP 


White Dwarf 


16, 


.9 


16, 


.6 


15. 


.4 


GD128 


Nov-Dec 05, Jan 06 


white 


OA, BA 


^Vhite Dwarf 


15 


.89 


15 


.82 






Hz2 


Nov-Dec 05 


V, 


b. 


u, uvwl, uvm2, uvw2 


OA, BA 


White Dwarf 


13, 


.86 


13, 


.81 






GD50 


Dec 05-June 06 


V 






PDU 


White Dwarf 


13, 


.98 


13, 


.82 






SA95-102 


March 05 


V, 


b. 


u 


EA, ZP, CT 


Landolt 


15, 


.622 


16, 


.623 


16. 


,785 


SA98-646 


March 05 


V, 


b. 


u 


EA, ZP 


Landolt 


15, 


.839 


16, 


.899 


18. 


,325 


SAlOl-278 


March 05, Nov 05 


V, 


b. 


u, white 


EA, ZP, CT 


Landolt 


15, 


.494 


16, 


.535 


17. 


,272 


SAlOl-13 


March 05, Nov 05 


V, 


b. 


u, white 


EA, ZP, CT 


Landolt 


15, 


.953 


16, 


.590 


16. 


,557 


SA104-244 


Feb-March 05 


V, 


b. 


u, white 


EA, ZP, CT 


Landolt 


16, 


.011 


16, 


.601 


16. 


.449 


SA104-338 


Feb-March 05 


V, 


b. 


u, white 


EA, ZP, CT 


Landolt 


16, 


.059 


16, 


.650 


16. 


.568 


SA104-367 


March- April 05 


V, 


b. 


u 


EA, ZP, CT 


Landolt 


15, 


.844 


16, 


.483 


16. 


,357 


SA104-443 


March- April 05 


V, 


b. 


u 


EA, ZP, CT 


Landolt 


15, 


.372 


16, 


.703 


17. 


,983 


SA104-457 


March-April 05 


V, 


b. 


u 


EA, ZP, CT 


Landolt 


16, 


.048 


16, 


.801 


17. 


,323 


PG1525-071b 


March 05 


V, 


b. 


u 


EA, ZP 


Landolt 


16, 


.403 


17, 


.133 


17. 


,268 


NOG 188 


Oct 05 


b 






GL 


Open cluster 














SA95 


July-October 05 


b 






GL 


Star field 














M67 


March 06 


b 






CL 


Open cluster 














SA104 


March 05 


b 






GL 


Star field 















Table 3. Optimum aperture results obtained from Hz2 (optical 
and UV filters) and GD128 (white filter). 



Filter 


Aperture radius (pixels) 


Range (pixels) 


V 


11.7 


2.6 


b 


10.8 


1.2 


u 


8.7 


1.8 


uvwl 


12.5 


0.9 


uvm2 


10.0 


0.5 


uvw2 


9.7 


0.3 


white 


9.8 


1.3 



radius for all filters including the white filter. The differences 
between the enclosed energy fraction within our adopted 10 
pixel radius aperture, and the optimum apertures we de- 
termined for the individual filters (Table [3]) are very small 
because this radius is already in the wings of the PSF where 
the count rates are low (Figure [5|. In the case of uvwl, 
which has the largest difference between the optimum aper- 
ture radius of 12.5 pixels, and our adopted aperture radius 
of 10 pixels, the difi'erence in the enclosed energy fraction is 
at most 2 per cent. 

6.2 Subtracting the background 

Two methods of determining the background count rate 
were considered: mean background, or clipped-mean back- 
ground. The mean background method averages the number 
of counts per pixel over the background aperture, whereas 
the clipped-mean background method excludes any back- 
ground pixels with count rates more than 3 sigma above 



the initial mean and then averages the counts over the re- 
maining background aperture pixels. The advantage of the 
clipped-mean background method is that it removes counts 
from any significant sources that may lie within the back- 
ground aperture. However, the disadvantage is that it pro- 
duces a consistently lower value because it removes the top 
of the distribution, which is assumed to be Gaussian. For 
very low rate backgrounds, as is often the case with UVOT, 
this assumption is not valid: most pixels contain either 1 or 
counts. 

Using the software package GAIA (starlink GAIA 
version 2.8-0) the clipped-mean method was consistently 
smaller than the mean method by 0.14 ± 0.03 ph/pix. An 
in-house IDL routine was also used, and also resulted in a 
consistently lower background value when using the clipped- 
mean method. This is acceptable when observations have 
high backgrounds, but will have a significant effect on count 
rates for observations with lower backgrounds and faint 
sources. Hence, a background limit of 10 ph/pix was set 
(where the measured background difference is 1.4 per cent of 
the total background), above which the clipped-mean back- 
ground method is applied, and below which the simple mean 
background method is used unless there is an obvious source 
that lies within the background region. 



6.3 Background region size 

For bright sources PSF wings can be seen that extend out 
to around 40 — 55 pixels. Investigation into the extent of the 
PSF wings was carried out by examining images of three 
white dwarf standard stars (Hz2, GD128 and WD112H-145) 
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1500 



1000 



500 




40 

Radius (pixels) 

Figure 5. Data for Hz2 in tlie b filter with a hardware win- 
dow size of 960 X 960. The dashed fine represents the number of 
counts in each consecutive ring for the source plus background 
measurements, the dotted line is the number of counts in each 
consecutive ring for the background only, and the solid black line 
is the number of counts in the rings once the background has 
been subtracted. The error bars on the background subtracted 
data are due to Poisson statistics. 



Table 4. Radius at which the count rate per pixel of the source 
becomes indistinguishable from the background. 



Filter 


Hz2 


GD128 


WD1121+145 




(pixels) 


(pixels) 


(pixels) 


V 


20 




15 


b 


50 




15 


u 


51 




21 


uvwl 


52 




45 


uvm2 


54 




52 


uvw2 


55 




48 


white 




30 


20 



of different magnitudes, with the three default hardware 
window sizes. 

For each case we measured the total number of source 
and background counts in 70 concentric rings around the 
source; each ring one pixel wide and varying in radius from 
1 — 70 pixels from the centre of the source. An example is 
shown in Figure[S] The background subtracted source counts 
are represented by the solid black line, and in this example 
the radius at which the total number of counts becomes in- 
distinguishable from the background, is at around 50 pixels 
radius. Table |4] shows the radius at which the count rate 
per pixel of the source becomes indistinguishable from the 
background. These results show that even in the case of the 
faint source, WD112H-145, very faint PSF wings are still 
sometimes observed. 

From these results the background region should be at 
least 55 pixels from the centre of the source, thereby avoiding 
any PSF wing photons in the background calculations. We 
used an annulus with inner radius 55 pixels from the centre 
of the source and an outer radius of 70 pixels to produce a 
large sample of background area with which to calculate the 
background contribution. 



6.4 Count rate calculations 

The method used to obtain corrected count rates for stan- 
dard star UVOT calibration observations is outlined in the 
following section. This is the procedure that was used for 
the calibration analysis, and is not necessarily a recipe for 
analysing GRB data. Please go to Section [14] for practical 
guidelines to analysing data. 

The raw image data from the UVOT was preprocessed 
using the standard Sm/t pipeline (UV0T2FITS v3.16). The 
first step of the pipeline is to construct a bad pixel map 
for each exposure and then to use this map to exclude pix- 
els from further analysis. Next, fk5 equatorial coordinates 
were applied to each exposure and the exposures were ro- 
tated so that north is up and east to the left. This coordi- 
nate system is applied using information from the Swift star 
trackers, which are accurate to five arcseconds. The pipeline 
corrects for this uncertainty by matching the star field in 
each exposure to sources in the HST Guide Star Catalogue. 
This aspect correction is accurate to 0.5 arcseconds (90 per- 
cent confidence radius) . The radial aperture size of 10 pixels 
meant that a fixed pattern correction (see Section (2]) did not 
need to be applied to these data. 

Following is a list of steps taken to determine the count 
rates of standard stars from the sky images produced by the 
Swtft pipeline: 

(i) We obtained an aspect-corrected image 
in SKY coordinates from the Swift archive 
(http:/ /swift. nasa. gov/docs/swift/archive ). 

(ii) We removed any exposures or observations that con- 
tained any exposure time anomaliefl. 

(iii) Where multiple exposures were taken in an observa- 
tion, count rates were calculated for the individual expo- 
sures, and then a weighted mean obtained. 

(iv) We obtained observed (source -I- background) counts 
using a 10 pixel (5 arcsec) radius aperture. 

(v) We obtained the dead-time corrected exposure time 
for each exposure to calculate the count rate (the EXPO- 
SURE keyword in the UVOT sky image files). 

(vi) We obtained the background count rate from an an- 
nulus with an inner radius of 55 pixels and an outer radius of 
70 pixels centered on the source. If the background level was 
below 10 ph/pix we used a mean background method, or if 
it was above this level we used a clipped-mean background 
method, as described in Section [6.21 

(vii) We corrected the (source -I- background) and back- 
ground count rates for coincidence loss using Equation [3] as 
described in Section [7] 

(viii) Finally, we subtracted the coincidence-loss- 
corrected background count rate from the coincidence-loss- 
corrected (source + background) count rate to obtain the 
coincidence-loss-corrected source count rate. 

Errors on the count rates were calculated using Pois- 
son statistics and were carried through the coincidence loss 
equations to produce final coincidence-loss-corrected count 
rate errors. 



Please see |http;//h eas arc. gsfc. nasa. gov/docs/swift/analysis/uvot_digest.html[ 
for more details on timing problems. 
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7 COINCIDENCE LOSS 

As described in Section (2] the UVOT detector suffers from 
coincidence losses at high count rate when multiple photons 
arrive at the same location on the detector during a sin- 
gle frame (see Figure [6} . The theoretical coincidence-loss- 
corrected count rate for a single-pixel detector is: 

^ _ -ln(l - aCra-wft) 

^theory p \^) 

where Ctheory is the incident count rate (in counts per sec- 
ond) and Cram is the raw observed count rate (also in counts 
per second) calculated using the dead-time corrected expo- 
sure time (keyword EXPOSURE in the UVOT image FITS 
file header), ft is the frame time (0.0110329s for a full frame; 
keyword FRAMTIME in the FITS file header) and a is the 
dead time correction factor (one minus the dead time frac- 
tion; 0.9842 for a full frame; keyword DEADC in the FITS 
file header). 

This theoretical coincidence-loss expression cannot be 
applied on a pixel-by-pixel basis to UVOT images because 
each count assigned to a UVOT pixel is determined by cen- 
troiding a photon splash over five physical CCD pixels (in a 
'cross-hair' sampling) . The effect of high coincidence loss is 
thus not only to lose counts, but also to reposition counts, 
since overlapping photon splashes within a single frame will 
be misplaced by the centroiding algorithm. 

We therefore introduce an empirical polynomial correc- 
tion to account for the differences between the observed and 
theoretical coincidence loss correction: 

/(x) — I + ai X X + 02 X + as X + a4, X x"^ (2) 

where x = Craw ft- Hence we have a full coincidence loss 
corrected incident count rate of: 

Ccorr — Ctheory f (^-^^ (3) 

where Ccorr and Ctheory are both in counts per second. The 
coefficients in Equation [5] were determined by a least-square 
fit to minimize the differences between the UVOT b magni- 
tudes of 361 isolated st ars in th e open cluster NGC 188, with 
the catalogue of Stetson et al. I ((2004) • The count rates were 
measured within 10 pixels as described in Section 16.41 The 
fit yielded the values ai = 0.066, 02 = -0.091, as = 0.029, 
and 04 = 0.031. Figure [7] compares the UVOT and Stetson 
photometry of the 361 stars in NGC188 after the polynomial 
correction has been applied to the UVOT data. This correc- 
tion was verified using UVOT observations of oth er photo- 
metric fields (SA95, M67, SA104) also studied by IStetson I 
(|2000l ). Note that the polynomial component changes the 
coincidence loss correction by less than 3 per cent for count 
rates up to 0.96 counts per frame (87 counts per second for 
a full frame), at which point the coincidence loss correction 
is a factor of 3.4. 

The coincidence loss correction here applies to isolated 
point sources and does not apply for crowded or extended 
sources, or when a smaller aperture is needed to maximize 
the signal to noise ratio. There is no general solution to these 
problems though we note that for faint sources, it may be 
preferable to maximize the signal to noise ratio with a small 
aperture, and accept some increased scatter in the small 
coincidence loss correction (see Section [2}. 

The finite number of frames in an exposure implies that 
the measured count rate follows a Binomial distribution. Be- 



Counts/Frametime 
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Figure 6. The size of the coincidence loss correction (the ratio of 
the incident to observed counts) is shown as a function of the ob- 
served dead-time corrected count rate for the nominal full-frame 
observing mode. The dotted line shows the expected relation for 
a single pixel device, while the solid line includes our empirical 
polynomial adjustment to the theoretical relation. The top axis 
shows the number of counts per frame. 
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Figure 7. The difference between UVOT and IStetson et al. I 
(200^) magnitudes for 361 isolated stars in the open cluster NGC 
188 is plotted against the number of counts per frame (bottom 
X-axis) and Stetson B magnitude (top x-axis). A polynomial cor- 
rection has been applied to the theoretical coincidence loss curve 
for UVOT to minimize any systematic trend with the Stetson 
magnitudes. For clarity, only selected Poisson error bars are in- 
cluded for counts per frametime of less than 0.3. 

cause of this, normal photometry packages which assume 
Poisson statistics will not give an accurate measure of the 
error, except at low count rates. Since the incident count 
rate must be derived by use of the non-linear coincidence- 
loss correction according to equation [TJ the error in the in- 
cident count rate, atheory, is given by (see iKuin fc Rosen I 
iiool) ): 

i 1 T / - I ^ raw ft-^ / .\ 

'^theory = -— r '■n(l±- ) (4) 

" aft [l — Craw ft) 
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Figure 8. Shaded plot showing how the sensitivity of the UVOT 
varies with position. The paler areas are more sensitive. The data 
are binned to make 4x4 bins with a minimum of 4 points per 
bin. 

where utheory is in counts per second, and araw ~ 
•\/ Craw{l — Crawft)/teiapsed (in counts per secoud) is the bi- 
nomial error in the measured count rate. For low count rates, 
the Poisson error is therefore a good approximation. For the 
highest incoming photon fluxes (more than 0.9 counts per 
frame) the upper error becomes larger than the lower error, 
but in most cases of interest, they are nearly equal. 

Based on the NGC 188 photometry, any additional sys- 
tematic error introduced by the coincidence loss correction 
is less than 0.01 magnitudes. This is confirmed in the com- 
parison data sets where no trend is seen with magnitude. 



8 POSITION-DEPENDENT UNIFORMITY 

The photometric sources were all measured near the middle 
of the detector, but any variation in sensitivity from position 
to position due to irregularities in the fibre bundle, photo- 
cathode etc., or any larger scale trend in sensitivity over the 
field of view clearly should be taken into account. 

Large scale sensitivity (LSS) variations over the area of 
the detector were measured using the count rate of a stan- 
dard star at a variety of positions on the detector. For this 
we used GD50, and several other stable stars to evenly sam- 
ple the whole field of view, giving 163 points in total. The 
standard deviation of the normalised count rate of all these 
points is 3.2 per cent, but for the central 1024 region only, 
the standard deviation is 2.0 per cent. A trend of count rate 
versus position could clearly be seen, so the data were binned 
spatially into 16 bins, with a minimum of 4 measurements 
per bin, to remove small scale variations. The LSS variation 
over the full field of view demonstrated by the binned data 
has a 2.2 per cent standard deviation with a peak to trough 
range of 8 per cent. Figure [8] is a shaded plot showing this 
variation. One corner has a higher sensitivity and the oppo- 
site corner lower. However, after binning, the central region 
of the detector (the central 1024 pixels, equivalent to the 
middle 4 boxes in the diagram) within which the photome- 
try measurements were taken appears to have more-or-less 
fiat sensitivity showing a standard deviation of 0.6 per cent. 

An upper limit to the small scale, or pixel to pixel sen- 
sitivity variation, is found from repeated measurements of 
the same object (GD50) at slightly different positions in the 
central area and is found to be at most 2.3 per cent standard 



Table 5. Ratio of observed to predicted count rates used to pro- 
duce an in-orbit correction curve. 



Filter 


Average count 


Standard 


Number of 




rate ratio 


Error 


Observations 


V 


0.698 


0.008 


13 


b 


0.725 


0.011 


12 


u 


0.753 


0.014 


11 


uvwl 


0.786 


0.009 


3 


uvm2 


0.740 


0.010 


3 


uvw2 


0.784 


0.009 


3 



deviation. This has been confirmed by the measurements of 
several stars used in the long term stability calibration, and 
also by the pixel-to-pixel variations in count rate in expo- 
sures taken with the on board LED. At present, it would be 
prudent to assume a 2.3 percent systematic error on indi- 
vidual UVOT photometry measurements. 

A more detailed analysis of large and small scale sen- 
sitivity and of the sensitivity changes over the lifetime of 
UVOT will be presented in a future paper (jBreeveld et al. I 
12007.) . 



9 IN-ORBIT EFFECTIVE AREA 
DETERMINATION 

The predicted effective area curves for each filter were cal- 
culated in units of cm^ by multiplying the UVOT reference 
instrument response (Figure [TJ by the ground based filter 
transmission curves (Figure [Sjl . This section describes how 
the curves were adjusted to produce the in-orbit effective 
area curves. 



9.1 In-orbit instrument response curve 

To obtain in-orbit effective area curves, the instrument re- 
sponse curve was adjusted using an in-orbit correction curve. 
The correction curve was produced using the standard star 
observations through the UV and optical filters. For each 
star, in each filter, the ratio of the observed count rate to 
the predicted count rate was computed. These ratios were 
averaged for all the stars in each filter. These average in- 
orbit/predicted ratios and their corresponding standard er- 
rors can be seen in Table [5] 

The correction curve was produced by fixing the ratio 
values at the central wavelengths of each filter (Table[TJ , and 
applying a spline fit. The correction curve was assumed to be 
constant longwards of the central wavelength of the v filter 
and shortwards of the central wavelength of the uvw2 fil- 
ter. Figure |9] shows the in-orbit correction curve. The shape 
of this curve, which represents the deviation from the ideal 
response, is not consistent with significant contamination of 
the optical surfaces by molecular material, which would have 
a greater effect in the UV than in the optical. This in-orbit 
correction curve was then multiplied by the reference in- 
strument response curve to produce the in-orbit instrument 
response curve shown as the solid black line in Figure [T] 
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Figure 9. In-orbit correction curve made by a spline fit to the 
ratios given in Table [5] The y-axis error bars show the standard 
error for each ratio; the x-axis error bars show the FWHM of each 
filter (Table [ij. 



9.2 In-orbit effective area curves 

The in-orbit effective area curve for each filter, except for 
white, was calculated by multiplying the in-orbit instrument 
response with the ground-based filter transmission curves. 
Figure [TO] shows the resulting in-orbit effective area curves. 

The white filter is a special case because it covers a 
very wide band, and so the shape of its effective area curve 
depends on the shape of the instrument response curve 
throughout the wavelength range. Unfortunately, it is also 
the one filter that does not have a well measured through- 
put, and so a different sequence was followed for the cali- 
bration of the white filter effective area. First, we combined 
the white filter predicted transmission curve with the in- 
orbit instrument response function to predict the count rates 
of the standard stars observed through this filter. We then 
compared these with the observed count rates, to find the 
average in-orbit correction factor for the white filter. This 
ratio was found to be 0.88 ± 0.05. Unlike the other optical 
and UV filters, the white filter in-orbit effective area curve 
was then calculated by combining the in-orbit instrument 
response curve with the theoretical white filter transmission 
curve and then multiplying by the extra factor of 0.88. 



10 PHOTOMETRIC ZERO POINTS 

The photometric zero point Zpt of each UVOT filter is de- 
fined as the magnitud e which is equivalent to 1 count per 
second l|Pogson Ill856l ). and is given by the equation: 



Zpt = nis 



~\- 2.5 log {Csource) , 



(5) 



where the magnitude of a source, and Csource 

is the count rate of an observed source after correction for 
coincidence loss and dead time. 

We have taken the o bserved spectrum of Vega from 
iBohUn fc Gilliland I (|2004l ') to define our UVOT magnitude 



system. Thus this spectrum of Vega represents mve 







in all filters. Transformations are provided in Section [TT] to 
convert from the UVOT system to the Johnson system. The 
zero points for the optical and UV filters were calculated 




Figure 10. In-orbit effective area curves for the UVOT filters. 



by standardising the count rates to the Vega spectrum. The 
expected count rate of each observed star {Cexp{i)) was cal- 
culated by folding its spectrum through the in-orbit filter 
effective areas (Section 19. 2|l . In the same way the spectrum 
of Vega was used to produce an expected Vega count rate 
{Cexp{vega)). The zero points {Zpt{i)) for each source in 
each filter were then calculated using: 



Zpt{i) = mvega + 2.5 log Cexp{vega) 



(6) 



The final zero point {Zpt) for each filter was calculated by 
averaging over all the observations in that filter. 

Figure [TT] shows the data used to produce the zero 
points for each of the UVOT filters. The final, mean zero 
point is shown with a dashed line, and the RMS scatter is 
illustrated with dotted lines to either side. The error bars 
shown on the individual points in the plots include the Pois- 
son error in the raw observed count rate, and the errors as- 
sociated with the stellar spectra used. The predicted errors 
on the Landolt stars in the optical filters were calculated 
using an estimate of the sy stematic error betwee n the Lan- 
dolt and Johnson system (|Menzies et al. I [ToU) added in 
quad rature with the Landolt colour term errors (jLandolt I 
Il992j ) . The predicted errors on the Landolt stars in the white 
filter were calculated to be 4 per cent due to the scatter 
of possible Pickles spectral matches and this was added in 
quad rature with the Landolt colour term errors (|Landolt I 
Il992i ). The Oke and white dwarf predicted errors in all fil- 
ters were calc ulated using a HST spectrophotometric error 
of 2 per cent (|Bohlin 20061 ) or an lUE error of 3 per cent 
(|Massa fc Fitzpatrick II2OO0I) . 

Table [6] shows the final zero points. For each filter the 
average error, which is the mean of the individual errors for 
each point, is given in Column 6. As expected, this average 
error is comparable to the RMS scatter about the mean, 
which is given in Column 5. The error shown in Column 4 is 
the standard error for each zero point, which is a measure of 
the error on the mean and is smaller than the RMS. How- 
ever, the standard errors for the UV zeropoints are based on 
only three data points in each filter, and the standard er- 
ror for the white zeropoint is based on only six data points, 
so these cannot be considered as valid estimates of the un- 
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Table 6. In-orbit zero points. The standard error (Column 4) is 
the error on the mean zero point for each filter; the RMS (Column 
5) gives the scatter; the average error (Column 6) includes uncer- 
tainties on individual measurements. Column 3 (Recommended 
Uncertainty) lists the recommended error for each zero point. 



Table 7. RMS error on residuals to colour fits. PS is Pickles star 
data fits and GRB is GR,B models data fits. Columns 5 and 6 give 
the minimum and maximum limits on the x-axis (of Figures 1121 
and I14I I respectively, which give the range of colours over which 
the transforms were calculated. 



Filter 


Zero Point 


Recommended 


Standard 


RMS 


Average 






Uncertainty 


Error 




Error 


V 


17.89 


0.013 


0.013 


0.04 


0.03 


b 


19.11 


0.016 


0.016 


0.05 


0.05 


u 


18.34 


0.020 


0.020 


0.06 


0.07 


uvwl 


17.49 


0.03 


0.013 


0.02 


0.03 


uvm2 


16.82 


0.03 


0.015 


0.02 


0.03 


uvw2 


17.35 


0.03 


0.012 


0.02 


0.03 


white 


20.29 


0.04 


0.023 


0.05 


0.04 



certainties on the zeropoints. For this reason, for the UV 
and white zero points, we recommend using the errors from 
column 6 in Table [6] as the uncertainties on the zeropoints. 
For the convenience of the reader we have listed the recom- 
mended zeropoint error for all filters in Column 3. 



11 COLOUR TRANSFORMATIONS 

Colour transformations are needed to convert from the 
UVOT magnitude system to any other system, in order to 
compare UVOT data with photometry from other sources. 
In this section we provide colour transformations to convert 
from the UVOT system to the Johnson system for a range 
of stellar spectra. We also provide transformations for a set 
of synthetic GRB spectra because Swift is primarily a GRB 
mission. 

The colour transforms from the UVOT ubv system to 
the Johnson U BV system w ere calculated for stars using 
Pickles spectra (jPickles II 19981 ). and for GRB model spectra. 
The model GRB spectra were generated assuming a power 
law continuum and varying degrees of dust extinction, with 
the form FJ^ exp"^ , where (3 is the GRB spectral index, and 
e — A\/ Ay is the relative exinction per unit wavelength. The 
wavelength dependence of the extinction was modelled on 
the SMC extinction law (Pei 1992), and GRB spectra were 
produced for /3 = 0, 1, 2, rest-frame visual extinctions = 
0.0, 0.2, 0.4, 0.8, 1.0 mag, and redshift in the range 0.3 < z < 
2.0. 

To produce the Johnson colours and magnitudes, the 
Pickles spectra and GRB models were folded through the 
Johnson response curves (described in Section [Sjl ■ To pro- 
duce the UVOT colours and magnitudes the same spectral 
models were folded through the UVOT in-orbit effective ar- 
eas to produce expected UVOT count rates [C source), then 
converted into magnitudes using: 



f'sourcc 



pt 



2.5 log(Csoitrce) 



(7) 



where maource is the magnitude of the source in the UVOT 
system, and Zpt is the zero point of each filter (listed in 
Table EJ. 

Figure [12] plots the difference between UVOT and 
Johnson colours against UVOT colours for the optical filter 
combinations. The stars in each plot represent the Pickles 
stars, and the triangles represent the GRB models. The 
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PS 


U-B 


u-b 


0.057 


-1.482 


1.871 


PS 


B-V 


b-v 


0.025 


-0.364 


1.935 


PS 


U-V 


U-V 


0.075 


-1.846 


3.558 


GRB 


U-B 


u-b 


0.078 


-1.380 


0.543 


GRB 


B-V 


b-v 


0.034 


-0.124 


1.483 


GRB 


U-V 


U-V 


0.102 


-1.505 


2.026 


PS 


B-b 


b-v 


0.020 


-0.364 


1.935 


-PS 


B-b 


u-b 


0.030 


-1.482 


1.871 


PS 


V-v 


b-v 


0.014 


-0.364 


1.935 


PS 


V-v 


U-V 


0.015 


-1.846 


3.558 


PS 


U-u 


u-b 


0.073 


-1.482 


1.871 


PS 


U-u 


U-V 


0.071 


-1.846 


3.558 


GRB 


B-b 


b-v 


0.001 


-0.124 


1.483 


GRB 


B-b 


u-b 


0.002 


-1.380 


0.543 


GRB 


V-v 


b-v 


0.004 


-0.124 


1.483 


GRB 


V-v 


U-V 


0.004 


-1.505 


2.026 


GRB 


U-u 


u-b 


0.011 


-1.380 


0.543 


GRB 


U-u 


U-V 


0.012 


-1.505 


2.026 



solid line in each plot shows the second order polynomial 
fit for Pickles stars; the dotted line shows the second order 
polynomial fit for the GRB models, and the dashed line 
shows where the colours in the two systems would be 
equal. The RMS error on the residuals, and the ranges for 
which the fits have been calculated, can be seen in Table [7| 
The colour terms in these Pickles star polynomial fits are: 

U-B=0.034[±0.007]+0.862[±0.007](u-b) 
-F0.055[±0.006](u-b)2 

B-V=-0.004[±0.004]+1.039[±0.011](b-v) 
-0.037[±0.007](b-v)2 

U-V=0.071[±0.010]-f0.899[±0.008](u-v) 
+0.018[±0.003](u-v)2 



The colour terms calculated from the second 
order polynomial fits for the GRB models are: 

U-B=0.086[±0.003]+0.886[±0.007](u-b) 
-l-0.050[±0.006](u-b)2 

B-V=-0.008[±0.001]+1.012[±0.003](b-v) 
-0.018[±0.002](b-v)2 

U-V=0.162[±0.002]-h0.904[±0.002](u-v) 
-h0.010[±0.002](u-v)2 

Again, the ranges of colours over which the transforms were 
calculated are given in Columns 5 and 6 of Table [7] 

The observations of one faint white dwarf star and nine 
Landolt stars were then compared with these transforms. 
Observed count rates (Cots) were obtained using the method 
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M 19.00 



1 6.85 - 



E 15.7.5 



1 6.70 - 



Star 



8,50 
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1 7.35 



Star 



Star 



Star 



Figure 11. Zero point calculations for the UVOT filters. Left column: v, b and u filters; right column: uvwl, uvm2 and uvw2 filters; 
bottom: white filter. Each point represents the zero point obtained from a single standard star, labelled individually on the x-axis. The 
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Figure 12. The differences between UVOT optical colours and 
Joiinson optical colours for Pickles stars (stars) and GRB models 
(triangles): UVOT u-b compared with Johnson U-B versus UVOT 
u-b (top plot); UVOT b-v compared with Johnson B-V versus 
UVOT b-v (middle plot); and UVOT u-v compared with Johnson 
U-V versus UVOT u-v (bottom plot). The solid lines represent the 
second order polynomial fits to the Pickles stars, the dotted lines 
represent the second order polynomial fit to the GRB models, 
and the dashed lines show where the colours in the two systems 
are equal. The polynomial fit parameters are given in Section [TTI 



described in Section [6.41 and converted into magnitudes us- 
ing Equation [7] Figure [TH] plots the Pickles star fits to the 
Johnson versus UVOT colours together with these observa- 
tional data. The error bars on the observed data show the 
Poisson error in the raw observed count rate on the x-axis, 
and the err ors associated with the Landolt colour terms in 
the y-axis l|Landoltlll993 ). This Figure shows that within 
the scatter of these observations, the fits produced with the 
Pickles stars agree with the observations. 

Figure [14] plots the difference between Johnson and 
UVOT magnitudes against UVOT colours. The stars in 
each plot represent the Pickles stars, whereas the triangles 
represent the GRB models. Figure [T4l shows that the Pickles 
spectra and GRB models follow different curves: the solid 
line in each plot shows the third order polynomial fit to the 
Pickles stars, and the dashed line shows the second order 
polynomial fit to the GRB models. The residuals to the 
fits are shown in the lower panel of each plot, and show 
good agreement within 0.05 magnitudes, apart from a few 
outliers. As in Figure 1121 the outliers are due to Pickles 




0.0 0.5 1-0 

oVOT b-v 




-2 -1 1 

,.,VOT u-v 



Figure 13. Johnson optical colours versus UVOT optical colours 
using observed data. The solid lines are the polynomial fit to the 
Pickles spectra shown in Figure 1121 and described in Section 1111 
The dashed lines show where the colours in the Johnson and 
UVOT systems are equal. 



spectra from stars with deep absorption features; we have 
included these in the plots to give the reader an idea of 
the appropriateness of the colour transformations for a 
range of star types. The RMS error on the residuals can 
be seen in Table [7] The upper and lower colour limits for 
each colour transform are given in Columns 5 and 6. The 
colour terms obtained from the Pickles polynomial fits are: 

B-b=0.021[±0.003]+0.005[±0.012](b-v) 

-0.014[±0.022](b-v)2-0.011[±0.010](b-v)^ 

B-b=0.011[±0.004]-0.011[±0.008](u-b) 

-0.008[±0.004](u-b)2-0.002[±0.004](u-b)^ 

V-v=0.029[±0.002]-0.009[±0.009](b-v) 

-0.037[±0.016](b-v)2-K0.017[±0.007](b-v)3 

V-v=0.026[±0.002]-0.014[±0.002](u-v) 

-0.005[±0.001](u-v)2-^0.002[±0.0005](u-v)^ 

U-u=0.042[±0.010]-0.130[±0.020](u-b) 

-l-0.053[±0.010](u-b)^-0.013[±0.010](u-b)^ 

U-u=0.069[±0.012]-0.093[±0.009](u-v) 
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+0.037[±0.007](u-v)2-0.007[±0.002](u-v)=' 

The colour terms obtained for the GRB model fits are: 

B-b=0.0i6[±0.0003]-0.009[±0.00f](b-v) 
-0.023[±0.00i](b-v)2 

B-b=-0.018[±0.0005]-0.045[±0.001](u-b) 
-0.014[±0.00i](u-b)^ 

V-v=0.023[±0.00i]-0.02i[±0.003](b-v) 
-0.005[±0.003](b-v)2 

V-v=0.0i0[±0.0007]-0.0i2[±0.0006](u-v) 
-0.0009[±0.0006](u-v)^ 

U-u=0.068[±0.003]-0.159[±0.007](u-b) 
+0.036[±0.006](u-b)2 

U-u=0.172[±0.002]-0.108[±0.002](u-v) 
+0.009[±0.002](u-v)2 
Figure [T5] plots the difference between the Johnson and 
UVOT magnitude versus UVOT colour for Pickles star fits 
in comparison with observed data. As before, the error bars 
on the observed data include the Poisson error in the raw 
observed count rat e, and the error s associated with the Lan- 
dolt colour terms l|Landolt Ill992l ). 



12 COUNT RATE TO FLUX CONVERSION 

To compute an accurate flux density it is necessary to fold 
the source spectrum through the effective area curves. How- 
ever, we have found for all but the white filter that there is 
not a strong dependency on stellar spectrum across a wide 
range of classes. Thus, in many cases an estimate of the flux 
can be obtained directly from the count rate and this can 
be particularly useful when plotting UVOT data with data 
from other instruments. Therefore a count rate to flux con- 
version factor has been calculated for each filter. Particular 
care should be taken where there may be significant absorp- 
tion or emission features in the wavelength range of the filter 
(Table m, and also with the conversion factor for the white 
filter, which depends strongly on spectral shape because of 
the wide wavelength range (1600-8000 A). 

The count rate to flux conversion for each fi lter was 
calcu lated, as described next, using Pickles spectra (|Pickles I 
Il998l ). and a subset of the GRB power law spectral models 
described in Section [TT] with redshifts ranging from 0.3 < 
z < 1.0. 

The effective wavelength (Ae f f) for each filter was cal - 
culated using the Vega spectrum (|Bohlin fc Gilliland Il2003 ) 
in the following weighted formula: 

where F„effa(A) is the Vega flux at a given wavelength. A, 
Earea{X) is the predicted effective area. The resultant ef- 
fective wavelengths are shown in Table |51 ft must be noted 
that tails of the UV filter transmission curves extend into 
the optical range, as shown in Figure (2] and thus the effec- 



Table 8. Effective wavelengths for each filter, for a vega-like spec- 
trum, as calculated using Equation |8] 



Filter 


Wavelength (A) 


V 


5402 


b 


4329 


u 


3501 


uvwl 


2634 


uvm2 


2231 


uvw2 


2030 


white 


3471 



Table 9. Average count rate to flux conversion factor results for 
the Pickles star spectra. The units on the conversion factor and 
errors are ergs cm~^ A~^. 



Filter 


Conversion 


RMS 




Minimum 


Maximum 




factor 






b-v (mag) 


b-v (mag) 


v 


2.61 X 10-1^ 


2.4 X 10- 


-18 


-0.36 


1.09 


b 


1.32 X lO"^'' 


9.2 X 10" 


-18 


-0.36 


1.09 


u 


1.5 X 10-16 


1.4 X 10" 


-17 


-0.36 


1.09 


uvwl 


4.3 X 10-1'^ 


2.1 X 10- 


-17 


-0.36 


0.1 


uvm2 


7.5 X 10-16 


1.1 X 10- 


16 


-0.36 


0.1 


uvw2 


6.0 X 10-16 


6.4 X 10" 


-17 


-0.36 


0.1 


white 


2.7 X 10-1^ 


7.9 X 10- 


-18 


-0.36 


1.09 



five wavelengths for the UV filters for very red spectra will 
be longer than those given in Table [S] 

A UVOT expected in-orbit count rate was calculated 
for each model spectrum in each filter by folding the spec- 
trum through the UVOT in-orbit effective area curves (Sec- 
tion I9.2|l . A flux value at the effective wavelength for each 
model spectrum in each filter was obtained by smoothing 
the spectrum into lOA bins to remove narrow spectral fea- 
tures, and then interpolating over four points around the 
effective wavelength. 

A count rate to flux conversion factor was then calcu- 
lated for each spectrum and averaged to produce a count 
rate to flux factor for each filter. The average count rate to 
flux conversion factor for the Pickles star spectra and GRB 
models can be seen in Tablesl9land ll0l respectivelv. The RMS 
error (a measure of the data scatter) on the average factor 
is also given in each table, along with the range of UVOT 
b-v colours over which the factors were calculated. Outside 
these ranges the conversion factors may not be applicable. 



Table 10. Average count rate to flux conversion factor results 
for GRB models. The units on the conversion factors and errors 
are ergs cm"^ s"! A~^. 



Filter 


Conversion 


RMS 




Minimum 


Maximum 




factor 






b-v 


b-v 


v 


2.614 X 10-16 


8.7 X 10" 


-19 


-0.12 


0.73 


b 


1.472 X 10-16 


5.7 X 10" 


-19 


-0.12 


0.73 


u 


1.63 X 10-16 


2.5 X 10" 


-18 


-0.12 


0.73 


uvwl 


4.00 X 10-16 


9.7 X 10" 


-18 


-0.12 


0.03 


uvm2 


8.50 X 10-16 


5.6 X 10- 


-18 


-0.12 


0.03 


uvw2 


6.2 X 10-16 


1.4 X 10- 


-17 


-0.12 


0.03 


white 


3.7 X 10-1^ 


4.9 X 10- 


-18 


-0.12 


0.73 
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Figure 14. The difference between Jolinson and UVOT magnitudes versus UVOT optical colours for Pickles spectra (stars) and GRB 
models (triangles): B-b versus b-v (top left) and u-b (top right); V-v versus b-v (middle left) and u-v (middle right); U-u versus u-b 
(bottom left) and u-v (bottom right). The solid lines represent the third order polynomial fits to the Pickles stars, and the dashed lines 
represent the second order polynomial fits to the GRB models. The polynomial parameters are given in Section ll II 
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Figure 15. The difference between Johnson and UVOT magnitudes versus UVOT optical colours for observed data. The solid line shows 
for comparison the polynomial fit to the Pickles spectra illustrated in Figure [HI and described in Section [TT] 
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Table 1 1 . Comparisons of U VOT measurements of optical mag- 
nitudes with literature magnitudes. The UVOT optical magni- 
tudes and colours have been transformed to the Johnson system 
using the colour transformations. 



magnitude 






mean 


rms 


no. 


no. 


UVOT-other 






difference 


error 


sources 


data 








(mag) 






sets 


^trans V 






0.019 


0.020 


81 


5 


btrans B 






-0.024 


0.016 


81 


5 


^trans U 






-0.036 


0.036 


65 


4 


{b - v)trans - 


(B 


-V) 


-0.028 


0.004 


59 


2 


{u - b)trans — 


{U 


-B) 


0.024 


0.007 


59 


2 



Figure 16. This plot shows how the white filter count rate to 
flux conversion factor varies with b-v UVOT colours. The stars 
represent the Pickles data, and the triangles represent the GRB 
model data. 



Figure 17. This plot shows how the white filter count rate to 
flux conversion factor varies with uvw2-v UVOT colours. The 
stars represent the Pickles data, and the triangles represent the 
GRB model data. 



The large error in the white filter factor is due to large 
diflterences in the convolution of blue and red spectra with 
the white filter wavelength range. The white filter factor for 
Pickles spectra and GRB models across the UVOT colour 
b — V and uvw2 — v are shown in Figures [16] and [IT] respec- 
tively, demonstrating the large scatter in the factor. In both 
figures the stars represent the Pickles stars and the triangles 
represent the GRB models. 



by uvotf lux) programs in another, while a test version of 
the CalDB incorporating the new calibration was employed. 

The tests directly compared measured optical magni- 
tudes with those published in the literature. The sources 
used for u, b and v were: two sets of observations of 23 stars 
in the GRB051021a field and 36 stars in the GRB051022 
field (from Henden ftp:/ /ftp. aasvo.org/public/grb ) ; 10 stars 
in the field of the SN2005am (as calibrated by iLi et al."l 
(120061')'): six sources in the GRB060218 field (calibrated by 
iHicken et al. I l|2006l )). In addition the photometry of stars 
in the field PG1 633-I-099B have been compared with mag- 
nitudes given bv lStetsonll|2000l ): this work is shown in Sec- 
tion [14T] 

Not all the measurements were made in the central 
region where the throughput is most uniform (see Sec- 
tion[H}. Those in the GRB060218 and SN2005am fields were 
distributed fairly randomly across the detector in individ- 
ual eposures. For the GRB051021a and GRB051022 fields, 
about half the stars fell outside the central 4 squares in 
Figure [8l but did not fall in the corners; the large scale sen- 
sitivity for the positions of the majority of these sources is 
estimated to vary from 0.99 to 1.02. 

The u, b, and v magnitudes were obtained from UVOT 
observations using the HEAsoft tools mentioned above. Af- 
ter converting these instrumental magnitudes into Johnson 
U, B and V magnitudes using the colour corrections (equa- 
tions in Section [TT|) . comparisons were made with published 
U, B, V photometry. 

In general, the differences in magnitudes as observed 
by UVOT and those obtained from references as described 
above are small (see Table lll|l . There is no evidence for a 
bias in any of the magnitude differences. The transformed 
UVOT colours (6 — v)trans and (u — b)trans also match well 
with the reference samples. 



13 VERIFICATION 

We compared UVOT optical photometry measurements 
with photometry from other instruments in order to assess 
the accuracy of the UVOT photometric system indepen- 
dently. Three different groups working independently, have 
performed the comparison using different data, and using 
software other than that used for the calibration work. The 
software used included the latest HEAsoft Swift uvotsource 
program in one case, and a combination of the uvotdetect 
and uvotmag (whose functionality has since been taken over 



13.1 Binned data 

The photometric calibration of the UVOT is based upon 
data taken in 1 x 1 binned (i.e. unbinned) modes. However, a 
large fraction of observations by the UVOT are made using 
2x2 binned modes, to reduce the telemetry volume gen- 
erated on board the spacecraft. "We have checked whether 
the binning mode has a significant effect on the photometric 
calibration, by comparing the count rates of a set of stars ob- 
served in 1 X 1 and 2x2 binned modes. The average difference 
in the count rates for 30 stars in the field of GRB060206, in 
the UVOT V band using a 5 arcsec radius aperture (5 pixels 
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Table 12. Aperture corrections are given here in terms of mag- 
nitude for apertures having radii from 2 to 4.5 arcsec, calculated 
with the uvotapercorr task Swift_Rel2.7(Bld21)_15Jun2007, us- 
ing the CalDB file swureef20041120vl02.fits. The 2.0 or 2.5 arcsec 
aperture radii are not generally recommended, for reasons given 
in Section |14I but may be the best option where the field is very 
crowded or the background is complex. 



Aperture 



(arcsec): 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


V 


-0.276 


-0.145 


-0.091 


-0.054 


-0.032 


-0.014 


b 


-0.327 


-0.176 


-0.111 


-0.065 


-0.037 


-0.015 


u 


-0.329 


-0.169 


-0.103 


-0.059 


-0.034 


-0.015 


uvwl 


-0.405 


-0.212 


-0.126 


-0.069 


-0.037 


-0.015 


uvm2 


-0.342 


-0.182 


-0.109 


-0.060 


-0.033 


-0.014 


uvw2 


-0.417 


-0.222 


-0.133 


-0.073 


-0.039 


-0.016 


white (b) 


-0.327 


-0.176 


-0.111 


-0.065 


-0.037 


-0.015 



this method is to be used, we caution against using an aper- 
ture smaller than 3 arcsec because the orbital variations in 
the PSF have a significant infiuence on photometry calcu- 
lated with smaller apertures. 

An alternative method that eliminates the problem of 
the small orbital variation, is to derive the aperture correc- 
tion on an exposure by exposure basis. This also automat- 
ically copes with the PSF filter dependency. Furthermore, 
during settling exposures, or very occasionally when there 
are attitude problems, the pointing can drift during the ex- 
posure, causing the images to be blurred. In these cases an 
exposure-specific correction is essential. The shape of the 
UVOT PSF also depends on the count rate of the source; 
bright sources (Js 10 counts s~^) have narrower PSFs than 
faint ones because of coincidence loss efi'ects. This count rate 
PSF variation could be compensated for by deriving sepa- 
rate aperture corrections for sources with different count 
rates. In practice both the orbital variation and count rate 
dependency changes the aperture correction by only a few 
percent, except for very high count rates, in which case a 
small aperture is not appropriate. 

We present here a suggested method for obtaining well 
calibrated data from UVOT, using a small aperture plus an 
exposure-specific aperture correction, and include an exam- 
ple of this method below. However, we do not wish to be 
overly prescriptive, or to tie the user to any particular soft- 
ware tools. The reader should not be afraid to experiment 
with different apertures or background regions for different 
situations. 

One method of aperture correcting for each exposure is 
as follows: 

(i) Identify a set of 5-15 isolated stars in the exposure. 
These stars should be bright enough that there is signal in 
the wings of the PSF, but not so bright that coincidence 
loss is distorting the PSF (i.e. between a few and about 10 
counts per second). There should be no neigbouring sources 
detectable within at least 10 arcsec to avoid using stars that 
are contaminated by light from other sources. 

(ii) Perform photometry on each of these stars with both 
the user-supplied aperture and the 5 arcsec standard photo- 
metric aperture. 

(iii) Subtract the magnitudes in the user-supplied aper- 
ture from the magnitudes in the 5 arcsec aperture to get the 
aperture corrections. 

(iv) The mean of the 5-15 aperture corrections thus ob- 
tained can then be used as the aperture correction for that 
exposure. To estimate the error in the aperture correction 
take the RMS of the individual values about the mean value. 

(v) The aperture correction is added to the magnitudes of 
the sources of interest measured in the user-supplied aper- 
ture in the same exposure. 

The error in the aperture correction factor for a single 
measurement is not simply the square root sum of the mea- 
surements at the chosen and the standard aperture. That is 
because the photons in the smaller aperture are counted in 
measurements with both apertures. It can be shown that the 
error in the aperture correction factor is the measurement 
error in the largest aperture. It is useful to use this error 
for deriving the weighted average from multiple measure- 
ments, especially if the stellar magnitudes of the reference 
stars vary. 



for the 2x2 binning) , is less than 1 per cent of the average 
count rate, which is less than the average measurement er- 
ror based on Poisson statistics. We conclude that there is no 
evidence for a significant difference in the count rates mea- 
sured in 1 X 1 and 2x2 binned modes. Thus this calibration 
should be relevant to 2 x 2 binned and unbinned data. 



14 PRACTICAL PHOTOMETRY 

The data analysis described earlier, in Section [6l concerned 
the method used to obtain an optimal calibration over a 
broad range of source strengths. If the reader were to fol- 
low this method (using a 5 arcsec aperture and a large 
source-free background annulus more than 27 arcsec from 
the source), they should obtain reasonable results. 

However, when doing photometry on point sources, par- 
ticularly faint ones, the size and shape of the extraction 
aperture needs to be optimized to take into account the size 
and shape of the PSF, the degree of crowding in the expo- 
sure, and the desired science. For isolated point sources the 
optimal aperture should maximize the signal-to-noise ratio. 
We have found that in most cases the maximum signal to 
noise is obtained with an aperture radius between 2.5 and 
4 arcsecs, but the actual value depends on the strength of 
the source, t he density of sou rces in the field and the back- 
ground level. iLi et al.l l|2006t ) similarly found an optimum 
aperture radius of 2.5 arcsec for unbinned UVOT data and 
3 arcsec for binned data for the fields they were working 
with. Their work focussed on comparing UVOT optical pho- 
tometry with standard fields, and their calibration is based 
directly on these apertures. 

Since the UVOT calibration is based on counts mea- 
sured within a 5 arcsec aperture, a correction must be made 
if a different aperture is used. The size of the correction 
depends on the PSF of the source of interest which in turn 
depends on the filter being used. The PSF is also observed to 
vary slightly throughout the orbit as the instrument changes 
temperature. The CalDB contains a set of 'average' PSFs for 
each filter which can be used to derive a correction; an 'av- 
erage' magnitude correction has been derived and included 
in Table [T^ for quick reference. The HEAsoft Swift software 
tools uvotsource and uvotapercorr use these PSFs to de- 
rive corrections for any specified user aperture. However, if 
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Figure 18. A comparison of magnitudes from UVOT data of the 
PG1633+099B field, obtained using the exposure-specific aper- 
ture correction method outlined in Section If 41 with those from 
IStetson I bOOdt) . 



14.1 Using an aperture correction 

The stars in the PG1633-I-099B field have been analysed 
and calibrated using the exposure-specific method outlined 
above. The results (shown in Figure llSp can be compared 
directly with Stetson's photometry giving confidence that 
the method works well. The mean offset between Stetson 
and UVOT photometry in this data set is < 0.01 mag, with 
a standard deviation of 0.06 magnitudes for stars brighter 
than 17th magnitude, and 0.18 magnitudes if all the stars 
are included. 

We tested the simple standard PSF approach as used 
by the uvotapercorr tool in comparison with the exposure- 
specific method described above. The comparison was done 
for two v-filter data sets. We used uvotsource (which calls 
uvotapercorr) to perform photometry on several isolated 
stars in each frame using both 3 arcsec and 5 arcsec cir- 
cular apertures. An aperture correction to convert from 3 
arcsec magnitudes to 5 arcsec magnitudes was computed for 
each frame and compared to the aperture correction derived 
from the v filter PSF given in the CalDB. We found that 
the aperture correction derived from the PSF for these two 
fields underestimated the brightness of the point sources by 
0.020 mag in the 00055751001 exposures and 0.019 mag in 
the 00276321001 exposures. These examples confirm the fact 
that the aperture corrections derived from the PSFs are re- 
liable to a precision of a few hundredths of a magnitude, but 
that aperture corrections specific to each exposure need to 
be computed if high-precision photometry is desired. Simi- 
lar offsets are found for different combinations of filter and 
aperture size. 



15 DISCUSSION 

The approach we have taken in this work is superior in a 
number of ways to that which was used to produce the ini- 
tial in -orbit UVOT photometric calibration (iBreeveld et al. I 
I2OO5I ). A full description of the earlier calibration is beyond 
the scope of this paper, but we note here three particularly 
significant improvements: 
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Figure 19. UVOT u-v compared with Johnson U-V for Pick- 
les stars (stars) as plotted in Figure [121 The solid line represents 
the second order polynomial fit to the Pickles stars given in Sec- 
tion 111! the dashed line prepresents the fit given in the initial 
cal ibration (CalDB released 20050805). The transformation given 
bv lLi et al. I l|2006h is shown by the dash-dot line. 



(i) The shape of the instrument response curve used in 
this work, which is based on measurements of the individual 
UVOT optical elements, is a much closer representation of 
the true instrument response than that used for the original 
calibration, which, although it was a direct measurement 
and consistent with the current instrument response, was 
not of high-precision, and was sparsely sampled in wave- 
length. 

(ii) The effective areas and zero points are determined for 
the same aperture for all filters (5 arcsec radius) , that is op- 
timised to minimize the enclosed-energy dependence on co- 
incidence loss and hence count rate. In our initial calibration 
the effective areas and zero points were based on different 
(and less-optimum) apertures for the UV and optical filters. 

(iii) The zero points for optical and UV filters are gener- 
ated using a single procedure in the new calibration, whereas 
in the initial calibration the optical zero points were ob- 
tained such that the mean difference was zero between 
UVOT and Johnson magnitudes of a group of standard 
stars: in effect calibrating the UVOT as though its opti- 
cal filters had a response identical to the Johnson system. 
Defining a UVOT-based magnitude system, and then facili- 
tating conversion to the Johnson system using colour trans- 
formations gives an inherently more accurate, stable, and 
understandable calibration. 

For comparison with previously-released colour transfor- 
mations the u-v to U-V fit for Pickles stars (as shown in 
the bottom panel of Figure I12|) has been plotted again in 
Figure 1191 but this time we have added the initial in-orbit 
calibration from CalDB 200508 05 (dashed line ) and also 
the c olour correction provided in lLi et al~l l|2006l ) (dash-dot 
line). iLi et al.~l (|2006l ) suggest no colour correction for the 
V and b filters, which would therefore be represented by 
the dashed horizontal line in the middle panel of Figure \n\ 
The GRB models and fit are not included in this compari- 
son plot because this is the first time colour transformations 
have been provided for GRBs. 
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For a full list of changes to the calibration ap- 
proach, procedures and observations, the reader is 
referred to the calibration documentation on the 
Swift website. The CalDB documentation and files 
containing the calibrations presented in this pa- 
per are available from the HEASARC CalDB at 
|http: / / swift.gsfc.nasa.gov/docs/heasarc /caldb /swift / 1 
in version 20070627 of the Swift/UVOTA CalDB. In the 
event that the reader wishes to compare new results with 
data processed with previous versions of the pipeline or 
calibration products, they are advised to reprocess the older 
data. The pipeline has been revised to take care of exposure 
time problems and to flag any remaining problems. 



15.1 Future photometric calibration 

Although we consider the photometric calibration of UVOT 
as presented here to be in good shape, there are a number of 
areas in which we are working to further improve and refine 
the calibration. For example, we plan to measure the sensi- 
tivity variation over the detector with a higher spatial reso- 
lution than was possible using the observations described in 
Section [51 For this purpose, we have made a sequence of 55 
observations of a dense star field (more than 300 stars) at 
a range of offsets and orientations, so that each individual 
star is observed at least 30 times, each at different locations 
on the detector. The high density of data across the detec- 
tor and large number of stars will allow us to model the 
detector response spatially to minimize photometric error. 
The goal is to use these observations to construct a flat-field 
calibration product which will enable us to achieve a pho- 
tometric response that is uniform to 1-2 per cent over the 
entire detector area. 

It would be desirable to increase the number of standard 
stars and the colour range on which the UV and white filters 
are calibrated. The main limitation is the lack of UV spec- 
trophotometric standard stars which are sufficiently faint for 
UVOT. This may be overcome in two ways: by using smaller 
hardware windows we can bring some brighter standards 
within the range of calibrated coincidence loss, and also we 
would like to make use of white dwarfs from the Sloan Digital 
Sky Survey l|lvezic et al. II2007I ) for which there are accurate 
models of the UV flux. 

We also plan to characterise the point spread function 
out to a large radius to enable, for example, reliable surface 
photometry. 

The detector response is expected to decline gradually 
with total radiation dose. From experience with XMM-OM 
we expect this to be at a level of 1-2 percent per year. Al- 
though this decline is not yet large enough to be measured, 
continued observations of standard stars will allow us to 
measure and calibrate the changing photometric response of 
UVOT. 

There is also more work to be done on various aspects 
of practical photometry, for instance, by developing an opti- 
mal sky-fitting algorithm to remove stars in the background 
region, given that the background is usually in the Poisson 
statistics regime rather than the normally- assumed Gaus- 
sian. Also, the parameterization of the PSF with count rate 
and with orbital variations will lead to the development of 
more accurate aperture corrections. 



15.2 Summary 

In this paper we have presented the in-orbit photometric cal- 
ibration of the UVOT and defined the UVOT photometric 
system. We have discussed factors which affect the accu- 
racy of the photometry. Any observational science depends 
critically on the calibration of the instruments; the UVOT 
calibration is itself based on trusted sources and data pro- 
cessed with the most up-to-date software. The comparison 
with other calibrated data sets gives us confidence that we 
have a robust, reliable and accurate system. 
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